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5— i ' 

By means of effective field theory techniques, we study the modifications of some properties of 

^A ' weakly coupled heavy quarkonium states propagating through a quark-gluon plasma at temperatures 

much smaller than the heavy quark mass, tuq. Two different cases are considered, corresponding 

l/~) . to two different hierarchies between the typical size of the bound state, r, the binding energy, E, 

the temperature, T, and the screening mass, mo- The first case corresponds to the hierarchy 
rng 3> 1/r 3> T ^> E 3> rriD, relevant for moderate temperatures, and the second one to the 
hierarchy tuq »T» 1/r ,mn 3> E, relevant for studying the dissociation mechanism. In the 
first case we determine the perturbative correction to the binding energy and to the decay width 
of states with arbitrary angular momentum, finding that the width is a decreasing function of the 
velocity. A different behavior characterizes the second kinematical case, being the width of s-wave 
D ■ states a non-monotonic function of the velocity, increasing at moderate velocities and decreasing 

£-*, ' in the ultra-relativistic limit. We obtain a simple analytical expression of the decay width for 

T ^> 1/r 2> mo 2> E at moderate velocities, and we derive the s-wave spectral function for the 
more general case T 3> 1/r , mo 2> E. A brief discussion of the possible experimental signatures as 
well as a comparison with the relevant lattice data are also presented. 
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• ■ I. INTRODUCTION 

o 

C*") ' Heavy quarks produced in the early stage of relativistic heavy-ion collisions are valuable probes of the medium 

that develops at later stages. They can be used to resolve its energy density, and eventually to understand which are 
the fundamental degrees of freedom of the system. Indeed, at sufficiently high energy densities matter should form a 
quark-gluon plasma (QGP) and the propagating heavy quarks should be capable to convey this information to us. In 
relativistic heavy-ion collisions (HIC) it is expected that the variation of the interaction between heavy quarks due to 
the creation of a hot medium should be observable. In particular, the Debye screening of the Coulomb-like potential 
between a heavy quark and a heavy antiquark was proposed in Ref. [l|, [2| as a dissociation mechanism, resulting in 
the suppression of the yields of heavy quarkonium, QQ, states in HIC. The low-lying heavy quarkonium states are 
considered as the most powerful probes because they are the only hadronic states that are able to survive above the 
deconfincment temperature. This is due to the fact that even at weak coupling, namely ignoring confinement, these 
states still exist. In addition, the vector states enjoy a rather clean dilcpton decay channel, which makes them easy 
to identify experimentally, although in HIC the corresponding background is not yet completely understood, see e.g. 

Suppression of charmonium states was first observed in Pb-Pb collisions by the NA50 Collaboration [4[ at relatively 
low center of mass energy (per nucleon), \/Snn = 158 GeV. However, in contrast with the naive Debye screening 
scenario, a complicated pattern emerged because of the various processes involving charm quarks, see e.g. [5[ for a 
recent experimental analysis. Charm quarks can indeed be produced not only by hard scattering at the early stage of 
the collision (prompt production), but also later on by collisions inside the QGP (non-prompt production) and their 
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density is sufficiently high that they can recombine in charmonium states. Bottomonium states give a much clearer 
signal, see e.g. the discussion in Q, because these states are very massive and can only be promptly produced by hard 
scatterings: the probability of generating these states within the QGP is extremely low. Moreover, since bottomonia 
are heavy and compact objects, they do not equilibrate with the expanding medium, and can be really considered as 
external probes. 

Recent experimental results by the CMS Collaboration [7H9| and by the STAR Collaboration [h| indicate a clear 
suppression of T states, meaning that the nuclear modification factor Raa, expressing the ratio of the yields of a 
state in HIC with respect to p-p collisions, is less than 1. The Raa decreases with increasing centrality and/or with 
increasing pt and is higher for the fundamental state, meaning that it is less suppressed. Indeed, the experimental 
results of Pb-Pb collisions at \JSnn = 2.76 TeV [9| seem to indicate sequential suppression of bottomonium states 
j^-flOf. in particular, integrating over centrality it has been found that Raa(Y(Is)) ~ 0.6, Raa(Y(2s)) ~ 0.1 and 
R AA (T{3s)) < 0.1. 

The effective field theory (EFT) techniques are very useful for the description of heavy quarkonia, because they 
are suited for handling systems with well separated energy scales. In the case of QQ in a thermal medium [lll - fl4| 
two distinct kinds of scales appear, namely the non-relativistic scales and the thermal scales. The non-relativistic 
scales are given by the mass of the heavy quark, itiq, the typical momentum transfer, 1/r oc rriQa s (a 8 = g 2 /tor is 
the QCD coupling constant), and the binding energy, E oc mqa 2 llal . We have assumed the weak coupling regime 
and identified the relative velocity between Q and Q with a s (see J16l . Il7j for reviews) . The relevant thermal scales 
to our analysis are the temperature T and the Debye mass tud oc gT. We shall discuss two possible hierarchies 
itiq ^> 1/r 2> T ^> E ^> mo and toq 2> T ~^> 1/r, mo S> E, which we shall refer to as Case I and Case II, 
respectively. If the bound state moves with respect to the medium, the EFT analysis becomes more complicated, 
because additional energy scales may appear [18| ■ We shall restrict our analysis here to the case of moderate velocities 
(»*1) for which no further scales are induced, so that Case I and Case II above can be safely addressed. However, at 
some instances we will push our results to the ultra-relativistic limit (v — > 1). Although this gives the correct results 
in the QED case [lg], and hence, we expect them to be sensible for QCD as well, one must keep in mind that they 
are on a less firm ground. 

Using EFT techniques it has been shown that, at least in perturbation theory, the dissociation of heavy quarkonia 
is not due to the Debye screening but to the appearance of an imaginary part in the potential [ll|, [12, [HI, [2(| ■ In 
other words, at high temperature heavy quarkonia disappear not because the binding energy vanishes, but because 
the thermal width becomes so large that the QQ state melts in the continuum. In QCD two different processes 
contribute to the thermal width: inelastic parton scattering — which is the dominant process for mo ^> E, and the 
gluo-dissociation process that corresponds to the decay of a color singlet state in a color octet induced by a thermal 
gluon; this process is dominant for mo -C E [2l|, [22| ■ The inelastic parton scattering is often referred in the literature 
as Landau damping, the reason is that this scattering is always mediated by an intermediate space-like gluon and can 
be related to the absorptive part of the gluon propagator. We shall also use this nomenclature from now on. 

The study of bound states propagating in the QGP at finite velocity is relevant for T states that are promptly 
produced in HIC and will cross the hot medium with a relative velocity v. In principle it might happen that heavy 
flavors are drifted by the expanding plasma. Indeed the PHENIX Collaboration J23l . |24| has observed a large v? 
of heavy-flavor electrons, suggesting that there is significant damping of heavy quarks while they travel across the 
medium. This picture has also received support from microscopic calculations of heavy quark diffusion in the quark- 
gluon plasma [25| . However, the elliptic flow of the T by the expanding medium should be negligible if the Debye 
length is larger than the typical distance between quarks, because heavy quarkonium at distances larger than its 
radius is colorless. Therefore, in both Case I and Case II the drift should be small and certainly less important 
for bottomonia than for lighter quarkonia, like the J/tp, which can be non-promptly produced and are expected to 
roughly comovc with the thermal bath. This is because before recombining both charm quarks have been drifted by 
the expanding QGP. 

In the first study of moving QQs performed in [26[, the hierarchy of scale of Case II was assumed, but only the 
real part of the potential was considered. The imaginary part of the potential was studied in QED in [181 ] . where the 
velocity dependence of the cylindrically symmetric real and imaginary parts of the potential were determined. In the 
present paper we extend the analysis of [18( in two directions. 

Regarding Case I, we consider QCD instead of QED; the main difference is that while in QED a proton and an 
electron will always form an electrically neutral state, in QCD a heavy quark and a heavy antiquark can be found in 
a singlet and an octet state, and this induces new terms in the computation. We determine the velocity dependence 
of the thermal width and of the energy shifts at the leading order. In particular, we find that at the leading order the 
energy shifts of the s-wavc states do not depend on the velocity (like in QED), but the energy shifts of all the other 
states depend on the velocity (unlike in QED). 

Regarding Case II, we extend the analysis of [18( by deriving an approximate analytical expression for the s-wave 
width as a function of the temperature and of v, valid for the particular hierarchy of scales T 3> 1/r ^s> mo 3> E. 



Moreover, considering the more general case, where T ^$> 1/r and mu 3> E but the product rmjj is arbitrary, we 
solve the corresponding Schrodinger equation numerically and determine the spectral representation of the two-point 
function. 

This paper is organized as follows. In Sec. [II] we discuss Case I, corresponding to the hierarchy tuq 3> 1/r ^> T 3> 
E ^> mo- We derive the expression of the width and of the energy shifts as a function of the temperature and of the 
velocity of the bound state. In Sec. lIIII we discuss Case II, corresponding to the hierarchy rriQ 3> T 3> 1/V ~ tud 3> E, 
this section is divided in two subsections, in the first one we do an analytical analysis of the case 1/r ^$> mr> while 
in Sec. IIIIBI we solve the Schrodinger equation numerically for 1/r ~ mj and determine the spectral function for 
various values of T and v. In Sec. IIVI we present a brief discussion of the observable consequences of the velocity 
dependent thermal width, we compare our results with existing lattice simulations and we draw our conclusions. In 
the Appendix [X] we discuss the framework used to take into account the effect of a moving thermal medium. In the 
Appendix IB"! we present some details and numerical checks of the procedure used in Sec. IIIIBI to derive the spectral 
amplitudes. 

II. CASE I 

The low-lying bottomonium states, T(ls) and rjb, produced at early times in relativistic heavy-ion collisions are 
likely to have a typical size, r, smaller than the inverse temperature during most of their evolution in the QGP. At 
intermediate times, the temperature is also likely to be larger than the binding energy. Having in mind this possibility 
we shall study in detail the particular case 

m Q > 1/r > T > E > m D • (1) 

This hierarchy of energy scales was considered in [13( for a thermal bath at rest. For a moving thermal bath, it was 
studied in full detail for the hydrogen atom in [l]|; here we generalize those results to QCD. 

The general formalism to deal with a moving thermal medium is reviewed in the Appendix[X] An important outcome 
is that in the bound state reference frame two additional energy scales should be considered: 



T + =T X —- and T_=TJ—-, (2) 

where v is the velocity of the medium with respect to the bound state. When v — >• 1 these scales are widely separated, 
a fact that must be taken into account in order to build the appropriate effective field theory (EFT) . For instance, in 
[18| the appropriate EFT was constructed for the case T + ^> 1/r ^g> T, which is different from the EFT obtained for 
the case T + ~ T ~ T_, valid for v <* 1. We shall mainly restrict ourselves to the latter case, and comment on v — >• 1 



results. Note, indeed, that the QED analysis of [18| shows that the results obtained with the EFT theory valid for 
T + ~ T ~ T_ and then naively extrapolated to the v —> 1 case, coincide with the ones obtained with the proper EFT 
with T + > 1/r > T (if no large log resummations are performed). Hence, our results may hold for the v — > 1 case as 
well. 

A. Matching between pNRQCD and pNRQCD HT L 

Since 1/r »Twe can take as the starting point the pNRQCD Lagrangian at T — [23, [2^] , which is obtained 
from QCD by sequentially integrating out energy scales of order tuq and of order 1/r, 



£ p nrqcd = flight + J d 3 r (tt {St [id -h s ]S + 0+ [iD Q - h ] 0} 



+ Tr{O t r- g ES + S t r-. 9 EO} + iTr{O t r- 5 EO + O t Or-.gE} + ... J , (3) 

where £i; g ht is the QCD Lagrangian for light quarks, g is the coupling constant, E is the chromo-electric field, S and 
O are the quark-antiquark singlet and octet fields respectively, and 

h a , o = ^+V s , + ... , F s = -^i , Vo= VA/2-C F )a 
rriQ r r 

(r = |r|) correspond to the singlet and octet Hamiltonians (the dots stand for I/vtiq corrections); hereafter Ca — 3 and 
Cp = 4/3. Thermal corrections to this Lagrangian are exponentially suppressed because the energy scales integrated 



out (m,Q and 1/r) are much larger than T. Note also that no dependence on the velocity appears at this stage because 
the velocity enters in the calculation through the scales T + and T_ in @ only. 

Because we are assuming that the binding energy, E, is much smaller than the temperature, we may integrate out 
energy scales of the order of T as well. If we do so, we obtain an EFT which is temperature and velocity dependent. 
This EFT was called pNRQCDhtl in [13J . where it was used in the case of vanishing velocity. We consider here the 
general case of non- vanishing velocity, following the analogous of the QED calculation developed in [18[ . At the order 
we are considering, the pNRQCDhtl Lagrangian is obtained from £ p nrqcd in Eq. ^ by replacing £u g ht — > £vHTL, 
and h a — > h a + 6V S , where £ v htl is the Hard Thermal Loop Lagrangian for a plasma moving with a velocity v [29| , 
and the potential SV S encodes thermal contributions to the singlet potential, which depend on the velocity as well. 
The expression of SV S can be obtained by a standard matching procedure, using dimensional regularization (DR) to 
regulate the IR divergences arising from the expansion T 3> E. In this case we have to consider the pNRQCD diagram 
in Fig. [TJ where the dipole vertices and the octet propagator can be read from ©, 



d D k 
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SV S = -ig 2 C F r l \ -jt—^j^ — L _ , . (fcgAj(fc) + fcifc,-Z>oo(fc))»j , (5) 



with Dfiv^k) the gluon propagator. Since T ^$> E and we use DR, the following expansion can be performed 

1 E-h , (E-h ) 2 



(E — h ) — fco + * e \ko — ie (feo — ie) 2 (ko — ie) 3 



(6) 



which corresponds to a temperature expansion, meaning that upon substituting this expression in Eq. ([5]) we can 
expand the thermal contribution of the singlet potential as follows 

8V S = 8V StT3 + 5V S . T 2 + SV s ,t + 0(a s r 2 E 3 ) , (7) 

where 5V s ,T n oc T n . In the Coulomb gauge 

SV S . T3 = -g 2 C F r 2 J £^k fa ^) 2,S(k 2 - ^)f B ^~^ , (8) 

6V S>T , = -g 2 C F n(E - h )r 3 J ^0 fa - **l) 2**(*g - 0/*( ^=j^ ) , (9) 

d ° k 1 ( \ MA ..,* fc 2w f \kp-v-k \ 



8V S , T = -g 2 C F n{E - hoYrj / -— 5 - ^ - 4/ 2^ - fe 2 )/ B ' L ' ' , (10) 



where Jb(x) is the Bose-Einstein distribution function, see the Appendix IA1 Notice that <$Vs 7^ vanishes because the 
integrand is an odd function of fc ; also 5V S} t vanishes because the integrand differs from an odd function by a S(k ), 
which forces the integral to be zero in DR. Hence, the only non- vanishing contribution is given by 

SV a>T * = -t^L ri{ E _ ho)rj (p>. + f{v)P P) , (11) 

where we use the following notation, already introduced in [18j . 

f(v) = 1 (t>(2 - v 2 ) - 2(1 - v 2 ) tanh- 1 («)) , (12) 

^(^). ^ = 5^-^). (13, 

We can manipulate ri(E ~ hg)rj in the same way as it was done in the v = case in [13| . obtaining a more compact 
expression 

2irC F a s T 2 TrN c C F a 2 T 2 r ( (r ■ v) 2 

3m 12 \ Htr 







FIG. 1: The singlet self-energy. The plain line represents the singlet field, the double line represents the octet and the curly 
line corresponds to a gluon. 

where N c is the number of colors. The correction to the singlet potential in the pNRQCDhtl Lagrangian is 
2irC F a s T 2 , nN c C F a 2 s T 2 r (^ , ^ , (r • vf M otf ^ , ^. , p3 , ,_g, 



W. = ^g + ° ^ s ^1 + /(«) + ^j-(l - 3/H) j + 0(a.r'J5 d , afr^) , (15) 

where the Q(a 2 r 2 T 3 ) contributions above arise from a s corrections to the diagram in Fig. [1] They have been 
calculated in [13| for the v = case. Differently from the hydrogen atom case [T^], the correction to the potential 
depends explicitly on the velocity. This could be expected from the fact that the Gromcs relation (which is deduced by 
assuming Poincare invariance) is violated at finite temperature [30(. However, as we shall detail in the next section, for 
the s-wave states the corresponding velocity dependence in the energy shifts cancels out at first order in perturbation 
theory. 

B. Computation in pNRQCDhtl and final results 

With the obtained pNRQCDhtl Lagrangian we can evaluate the thermal corrections to the binding energy and 
to the decay width of the various hydrogen-like states. Since we are using perturbation theory, we shall assume that 
the wave-functions are given by the unperturbed hydrogen atom solutions, which can be identified by the principal 
quantum number, n, the angular momentum, I, and the magnetic quantum number, m. For a given heavy quarkonium 
state the binding energy at the leading order in the perturbative expansion is given by 

E nim = K + (nlm\mV s + $Y, s \nlm) , (16) 

where E^ is the binding energy taking into account only the Coulombic part of the potential and S s is the self-energy 
of the singlet component of the heavy quarkonium. Clearly, the Coulombic part of the potential does not contribute to 
the decay width, which is nonzero only because of the thermal corrections, and at the leading order in the perturbative 
expansion it is given by 

Tnim = —2(nlm\$s5V s + QTi s \nlm) . (17) 

The singlet self-energy can be determined computing the diagram in Fig. Q] but this time in pNRQCDhtl- In order to 
properly take into account the moving thermal bath, the boosted Bose-Einstein distribution function has to be used, 
and since T ^> E, we expand 
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v 
hence the self-energy can be written as follows 



\/l -v 2 J \k -v-k\ 



■ 2~ fv uM / d D k TVT^5(k 2 -k 2 ) i ( kikA , _, 2cn3 , Mn . 

S s = -ig z C F n{E - ho) 2 / tttwt n *-Ti TB T~\ i — T~ % ~ IT r J + 0{a s r 2 E A ) . 19 

J (2ir) D \k -v-k\ (E-h o )-k + ie\ k 2 J 

This integral is very similar to the one evaluated in the QED case [18| and can be written as follows 

S s = -ig 2 C F n(E - hof^JijTj , (20) 

where 9£Jy is given in Eq. (48) of (l8j ] . 



By manipulating ri(E — h ) 2 rj we obtain the following result 
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+ 0(a s r 2 E 3 ), (22) 



and at the considered order it is pure imaginary. On the other hand, the correction to the singlet potential in Eq. (fT5 
is real, with no imaginary part. Therefore, at the leading order 
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The correction to the binding energy in the regime 1/r > T > £ > mp is given by 
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where (l"l'm"m'\lm) are the Clebsch-Gordan coefficients and 
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(23) 



0{a s r 2 E 3 ,a 2 r 2 T 3 ), (24) 



(25) 



is the expectation value of the radial position operator in the hydrogen atom, with clq = l/mQCFOt s the Bohr radius. 
In the expression of the binding energy we can distinguish three different contributions, corresponding to the three 
terms in the square bracket of Eq. (|24[) . The first one is an overall energy shift, independent of the quantum state. The 
second term is a shift of the binding energy that removes the degeneracy in I associated to the so-called "accidental" 
symmetry of the hydrogen atom. The third term in the square bracket is an energy shift that depends not only on n and 
I, but also on m and it is thereby related to the breaking of rotational invariancc. There exists a privileged direction 
corresponding to v, thus the binding energy depends on the relative orientation between the angular momentum and 
the velocity. For s-wave states, having zero angular momentum, there is not such a dependence and this contribution 
to the binding energy vanishes, 
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(26) 



and as already anticipated it does not depend on the velocity of the plasma. The latter result is surprising, because 
one would have naively expected that a v dependence should arise because for the moving bound state the effective 
temperature depends on v, see the Appendix [X] However, our calculation shows that this is not the case, at least at 
the first order in perturbation theory 

Regarding the width, from the expression of the self energy in Eq. (|22[) , we obtain 
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is a negative and decreasing function of v. It can be easily shown, using the expression above, that for any state the 
width is a decreasing function of the velocity, vanishing for v — > 1, meaning that an ultra- rclativistic velocity has the 
effect of stabilizing the system. This behavior is due to the Vl — v 2 prefactor in Eq. (|2"7| , which can be traced back 
to the expansion of the boosted Bose-Einstein distribution function in Eq. (fl~8|) and is therefore due to the "Doppler 
shift" of the temperature, see the Appendix [S] As we shall see in the next section, an analogous behavior is obtained 
for the hierarchy of energy scales of Case II in the ultra-relativistic limit, although the microscopic description appears 
to be different. 



In the expression of the width we can further distinguish two contributions, corresponding to the two terms in the 
square bracket in Eq. (|27p . Both are velocity dependent, but the first one does only depend on n, meaning that it 
originates from terms that do not break the rotational and the accidental symmetries. The second term depends on 
all the quantum numbers and vanishes for s-wave states. Thus, for s-wave states the width simplifies to 

F™ = ^ T f^ (_*K + _J^_ + «^T) log (1±A + (ay E z,ayT*) , (29) 

6v \ rriQ mQa n z 8 / \1 — V J 

which, as observed above for the general case, is a decreasing function of the velocity, vanishing for v — ¥ 1. 

III. CASE II 

The dissociation of heavy quarkonium is expected to occur for 

T>l/r,m £ )> J E, (30) 

as it happens in a thermal bath at rest. In the color-screening model introduced in Jl| the dissociation takes place 
because the number of bound states supported by a Yukawa potential decreases with the range of the potential, 
which is proportional to the screening length (l/m£>). The effect of screening becomes important when the screening 
leng th is of the order of the size of the system (1 /mp ~ r) . However in the actual real-time potential computed in 
fl9| (and confirmed by the EFT computations [IJ, UJ^) the dissociation takes place because the potential develops 
an imaginary part (Landau damping), and bound states turn into wide resonances as the temperature increases. 
This effect becomes important at a parametrically different scale, l/(Tm| 5 ) 1 / 3 ~ r (up to logarithms) \U\ . It is 
then particularly interesting to address the question whether the mechanism of dissociation (screening versus Landau 
damping) remains the same when the bound state moves with respect to the thermal bath. 

An EFT study of this situation in the QED case for muonic hydrogen with massless electrons was already performed 
in [l8[ (see [Tj] for the v = case) . For heavy quarkonium the results are analogous and can be read off that study by 
changing the value of the Debye mass to the QCD case and by correcting for trivial color factors, as it was already 
pointed out in [l8( ■ We briefly review the basic steps of the derivation below. 

• Since mq 3> T, the starting point of the calculation can be the NRQCD Lagrangian at zero temperature [la ]. 
By integrating out the temperature scale we arrive at the NRQCDhtl, an EFT whose Lagrangian is the sum of 
NRQCD for the heavy quark sector (with thermal, velocity independent, corrections to the heavy quark mass) 
[111 ], and the HTL Lagrangian for gluons and light quarks, which now depends on the relative velocity v between 
the thermal bath and the bound state [29j . 

• We can also integrate out the scales 1/r and mc, which leads to pNRQCDhtl- Since gluons and light quarks 
in the HTL Lagrangian develop a mass gap of the order of mc, this effective theory does not contain them as 
explicit degrees of freedom, and hence it reduces to a singlet heavy quark-antiquark field interacting through a 
potential V s (which depends on v as well) . The main conceptual difference with respect to the case discussed in 
the previous section is that now, in general, the thermal contributions cannot be considered as a perturbation 
in the potential. 

The potential V s coincides with the one that was computed numerically in Sec.V of [l8| . In that paper qualitative 
arguments were put forward on how the dissociation mechanism is modified when the velocity of the bound state with 
respect to the plasma increases. We will quantify those arguments here, by focusing on the effects of this potential 
on the physics of the Is state. We shall discuss two different cases. 

1. We consider the particular case 1/r ^$> mp. The thermal contributions can still be considered as a perturbation 
to the Coulomb potential. This allows us to compute the leading thermal corrections to the decay width 
(almost) analytically and to derive some explicit expressions for the velocity dependence. In particular we can 
parametrically estimate how the dissociation temperature depends on the velocity if 1 — v 3> mpaQ. Then, alike 
the v = case, the dissociation mechanism is dominated by Landau damping effects. This has to be contrasted 
to the Case I in Sec. [II] where the decay width is entirely due to gluo-dissociation. 

2. We consider the general case in which the relative size between 1/r and mp is left arbitrary (1/?' ~ tud) 
and compute the spectral function. Although the concept of dissociation temperature is useful for qualitative 
estimates, there is no universal definition for it, and hence it is of limited usefulness for a quantitative comparison 
of our results with other approaches. On the contrary, the spectral function is a well defined quantity so that 



our results can be straightforwardly compared with those obtained by different methods, in particular by lattice 
computations. Furthermore, it is directly related to a physical observable, the thermal dilepton production rate 
[3l|,|32|. In current HIC experiments, however, the heavy quarkonium states are not expected to be thermalized, 
but rather to act as hard probes of the medium. In the spectral function a bound state with zero decay width 
appears as a delta function whereas scattering states produce a smooth curve. The spectral function allows us 
to observe all the intermediate situations, which manifest when changing the thermal bath temperature and 
velocity. 

For the expressions of the coupling constant and the Debye mass, we will use the following parameterization (we 
set N c = 3,N f = 3), 

a s = a s (l/a ) m D = ^ = = J 4xa s {2-KT)T*{N c + N f /2) , (31) 

v /31og(2 7 rT/A M g) V 

where the MS renormalization scheme has been adopted with A^jg = 250 MeV; we also fix mq = 4.881 GeV and the 
Bohr radius of T(ls) is given by ao — 0.74 GeV -1 , both values are taken from [33|]. This choice is motivated from the 
fact that computing higher order corrections to the potential would introduce a dependence on the renormalization 
scale of the type log™(r/i); on the other hand, computing higher order corrections to the Debye mass would introduce 
terms proportional to log™ (2ttT / fj,) . 

A. The particular case T >• 1/r 3> run » E 

In this case the potential can be considered as the Coulomb potential plus a perturbation, and hence the following 
formula provides a good approximation to the decay width of a s-wave state 

T s- W ave = _ 2 <„00|9V.(r)|n00) = - J d 3 r|^„(r)| 2 J -^(e lkr - l)A s (k,v) , (32) 

where ■0„(r > ) is the wave- function for a s-wave state in the Coulomb potential and Ag is the symmetric part of the 
00 component of the gluon field propagator in the Coulomb gauge, which has been computed in [l8| 

A (i. „\ 87r 2 a s C F Tmlf(v,e) 

As{k > v) = k(ki + m i g ( z , vmk i + m w( z , v ))> (33) 

and depends on v, k = \k\ and on 9, the angle between the vectors k and v. In the above equation 

2(1 — t^sm 0)°l 2 
and we have made the dependence on the Debye mass explicit by defining g(z,v) = IIr^, v)/rrr D , with 

(35) 



\J\ - v" 1 



sm 



and Ur(z,v) is the retarded self-energy of the 00 component of the gluon field in the Coulomb gauge with kg = 0, 
which was first computed in [26I] 1 . In the reference frame where the bound state is at rest 

U R (z,v) = a(z) + z p^, (36) 

1 — v z 

where 

' z + 1 + ie 



a ( z ) = — 



^-(^-l)-ln 



2 V z - 1 + ie 



(37) 



1 In l26t there is a misprint in the first line of Eq. (8), in which the global sign must be the opposite. 



and 

'z + 1 



b(z) = & 1) 



a{z)-m D (l-z') l--bi 



2 Vz-1 



(38) 



In principle we can obtain the decay width for any s-wavc state, but for illustrative purposes we shall focus on the 
ground state (n = 1). It is convenient to start the computation in Eq. (|32[) by first performing the integration over r, 

T s-u,aue = 2a s C F Tm 2 D f dcos6f(v,6) f°° —, *—. -^ ^— -- ( 1 J5-5— | > ( 39 ) 

1 D i_i ' 'Jo (k 2 +m 2 D g(z 1 v))(^+m 2 D g*(z,v))\ ( 1 + ^i) 2y /' 

where we have switched to cylindrical coordinates in momentum space and performed as well the trivial integration 
over the azimuth angle. 

The above expression can be numerically integrated, but we first obtain an approximate expression valid at moderate 
velocities. In this case, for any angle, g(z, v) is of order 1 and there are only two scales in the previous integral, 1/ao 
and mo- Moreover they fulfill the relation 1/ao ^S> tud, so that the technique of threshold expansion [34[ can be used 
to work out the integral, thus obtaining 



T s^wave = as c F Tm 2 D a 2 f d cos 0/ CM) (log 



( 2 \ 1 g(z,v)log(g(z,v))-g*{z,v)log(g*(z,v)) 



\m D a J 4 2(g(z,v)—g*(z,v)) 

+ O((m D a ) 2 )) . (40) 

This equation can be further simplified by taking into account that the log I m 2 a J is logarithmically bigger than the 
rest of the terms in the parenthesis. With this approximation we arrive at the following result, 

r s-wave 2a s C F Tm 2 D al / 2 \ 
1 a/I - 1? \m D a J ' l ' 



or, equivalently, 

T{- wave (v) 



r s- W ave^ = y ^Y 



(42) 



which holds up to 0(l/log(m_Dao)) accuracy and is independent of the heavy quark mass and of the temperature. 
Then, in the regime T 3> 1/r ^$> mu, the decay width increases with the velocity, as far as it remains moderate 
(v & 1). Note that this behavior is opposite to the one observed in the regime 1/r ^> T ^> E ^> mo in Sec. Ill Bt sec 
Eq. (|29|) . If we take into account that m 2 D is proportional to T 2 , the decay width at temperature T and velocity v is 
the same as the one that we would observe at v = but with 

r ^=(i— !^w< (43) 

provided that T v ~ T 3> 1/r 3> run- 

Results beyond the logarithmic accuracy of ([32)1 can be obtained by evaluating (|40)) . or even better (|39|) . The 
numerical values of T(v)/r(Q) for the T(ls) state are reported in the left panel of Fig. [21 for three different temper- 
atures, together with the approximate expression (|42[) . solid black line, for < v < 1. The approximate expression 
correctly reproduces the numerical values for v < 0.5, but for larger values of v the ratio of the width decreases and 
becomes temperature dependent, departing from Eq. (|42")) . This is due to the fact that for v — > 1 further scales must 
be considered (T + ^> T ^> T_) and ([4^)1 does not hold. This expression relies on the fact that 1/a 2 , ^> m D \g(z,v)\, 
which not always holds for v — > 1, even if 1/a 2 , 3> m|j does. 

In order to ascertain the reliability of the expression in ([3T))) and the origin of the difference between (|39[) and (|42l) , 
let us scrutinize the velocity and angular dependence of m 2 D g(z,v). We can interpret the square root of m 2 D g(z,v) 
with positive real part [la. l26l| as the velocity dependent Debye mass itid{v,9) (itid(v,9) should not be mistaken by 
the parameter mu that we used before, they coincide at v = only and have the same size for moderate velocities 
m * 1 only). In Fig. [3] we report the plots of #t[mr,(v,9)]/mD and of $s[mr)(v,0)]/mr) as a function of 9 for 
v = 0.1,0.5,0.9,0.99. The real part is peaked at 9 = n/2, corresponding to a vanishing value of the imaginary part, 
which is instead peaked at a value of 9 that with increasing v approaches it/2. For 9 no 7r/2, mr,{v, 9) is small and the 
imaginary part is of 0(m£>) for any value of v, meaning that the bound state can be approximated with a Coulombic 
wave function, and therefore in this region (|32l) represents a good approximation. Moreover, for v < 0.5 both the 
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real and the imaginary part of m 2 D g(z 1 v) arc of 0(itid), irrespective of the value of (9, and therefore the approximate 
expression (J32) is reliable. This approximation is still qualitatively good up to v ~ 0.9, although the increased value 
of m,D{v, 8) for 8 ~ 7r/2, suggests that the quantitative agreement might be lost, as indeed can be observed in the left 
panel of Fig. [2] 

The angular region that may jeopardize the perturbative expansion about the Coulombic wave function is only 
for 8 ~ ir/2 and for v > 0.9. Indeed for 9 = tt/2 the real part of the Debye mass has a peak, and for 8 ~ tt/2 
the imaginary part is large. However, this angular region is small. In order to quantify this region we consider 
Sj[m|j(u, 6)] — 2^R[mi){v, 8)]^s[itid(v, #)], which has two maxima for 



arccos 



1 



t;r- V 11 + v 2 - \/oV35 + I0v 2 + 3v 4 
2v 



(44) 



and therefore for v > 0.9 the angular region around 8 = tt/2 where the real and the imaginary parts are large is given 
by 



A6 = \d+ - 6-\ ^ y/2(l-v). 



(45) 



which clearly shrinks to zero for v — > 1. In order to clarify that the contribution of this angular region is small, we 
plot in the right panel of Fig. [2]the thermal correction to Coulomb binding energy (100|3?(5V A (r)|100), SV being the 
thermal contribution to the singlet potential, normalized to the Coulomb binding energy. This quantity should be 
small for (f39|) to be reliable, as it turns out to be the case. 

Note, however, that the angular region A8 gives for v — >• 1 the largest contribution to (|39p . and in this case 
g{z,v) ~ b(z)/(l — v 2 ). Since m 2 - ) \g(z,v)\ S> 1/ajj, the approximate expression in (|1(J|) does not hold anymore, in 
agreement with the results in the left panel of Fig. [21 We find in this case, using the same techniques of integration 
by regions, that for v — > 1 the decay width goes to zero like a s T\Jl — v 2 whereas the energy shift goes to a constant, 
consistent with the results reported in the right panel of Fig. [2] 

The velocity dependence of the width in the ultra-relativistic limit is similar to the one discussed in the Case I, in 
Sec lIIB! although the microscopic mechanism is different, Landau damping in the present case, and gluo-dissociation 
in that case. The reason for the decreasing of the decay width is probably related to the fact that a moving bound 
state feels a plasma with a non-isotropic effective temperature 



T oS (8,v) 



tVT 



1 



(46) 



see the Appendix \K\ for more details. Actually, the effective temperature is higher than T in the forward direction, 
and lower than T in the backward direction, thus it is not obvious that the width of the moving bound state should 
increase — or be modified at all — when the bound state moves with respect to the thermal medium. However, in 
the ultra-relativistic case the effective temperature is almost everywhere less than T, see Fig. [71 and it is higher than 
T only in a narrow region 6 ~ 0. According to the previous discussion this angular region does not give the leading 
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FIG. 2: (Color online) Left panel: Width I»/T(0) from Eq. |@2) (solid line) and I»/T(0) of the T(ls) from Eq. ([39]) at 
a few values of temperature, as a function of the velocity v. Right panel: Expectation value of the real part of SV/Ef as a 
function of the velocity, SV being the thermal contribution to the singlet potential and E\ the Coulomb binding energy for the 
Is state. 
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FIG. 3: (Color online) Real (solid red) and imaginary (dashed blue) parts of the velocity dependent Debye mass as a function 
of the angle between the vectors k and v, for v = 0.1, 0.5, 0.9, 0.99. 

contribution to the width, which is instead dominated by the 8 ~ n/2 region, where the heat bath is effectively cold. 
Thus, a velocity close to 1 tends to stabilize the system. 

The presence of an imaginary Debye mass for any v > can be related to the collisionless transfer of energy between 
the heavy quarks and the gauge fields. An accurate description of this phenomenon would require the discussion of the 
propagating modes, but an imaginary part of the Debye mass does in any case signal an instability. This phenomenon 
is akin to the plasma instabilities generated by a charged current in a plasma. A similar result was indeed obtained 
in [35|, [36|] where the destabilizing effect of a single heavy quark propagating in a thermalized QGP was studied. 



B. 



The general case T 3> 1/r , mo 2> E: the spectral function 



The s-wave spectral function was computed in [33, 38] for the case at rest (v = 0). 

The procedure developed in Ref. [37| can be easily generalized to a moving bound state. We shall use the expression 
of the potential determined in [18| . which takes into account the relative velocity, v, between the bound state and the 
expanding plasma, and consider that the system has cylindrical symmetry, with symmetry axis the direction of v. 

The formalism introduced in [37J can be generalized to cylindrical coordinates, resulting in the following expression 
of the spectral function 



p(uj) = lim / dt {cos(wt) 5R[^(i, r, z)} - sin(wi) ^[ip(t, r, z)]} , 

r— >0 

2->0' 



n 



where ip(t, r, z) = u(t, r, z)/r, and u(t, r, z) is the solution of the Schrodinger equation 

i a 2 i a 2 i 19 



id, 



d 2 i 



d 2 



rriQ o*z ttiq o*r rriQr^ rriQr or 



lh 



V(r,z))u{t,r,z) = 0, 



(47) 



(48) 
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with the initial condition u(0, r, z) = —6N c rS 2 (r)8(z). 

In order to numerically handle the two-dimensional Schrodinger equation, we use an operator- splitting method, 
namely we split the differential equation in two differential equations each containing derivatives with respect to only 
one variable (z or r). Therefore, the whole Hamiltonian is divided in two pieces, H = Hi + H 2 , where 



Hi = 
H 2 = 



[^Wz- V ™ 



m,Q d 2 i 



1 



1 d 
rriQr dr 



(49) 
(50) 



and we then solve the corresponding Schrodinger equation recursively in a discrete space-time, see the Appendix [E] 
for more details and for a check of the numerical code. 

In Fig. [4] we report the spectral functions at vanishing velocity for certain values of the temperature as a function 
of lu/ttiq, where u> <C tuq is the non-relativistic energy (w = corresponds to a relativistic energy of Imq). For the 
sake of comparison, we also report (right panel) the results we obtain for the choice of a s and mn of Refs. [371 139| . 
which are in quantitative agreement with those obtained in j37|. At T — 250 MeV the spectral function is given by 
the superposition of a peak, corresponding to the T(ls) bound state and of a continuum. The bound state has a 
thermal width which is determined by the imaginary part of the potential and is dominated by the Landau damping. 
The width of the bound state increases with increasing temperature, and correspondingly, the contribution of the 
continuum increases. At T ~ 500 MeV no peak of the spectral function is visible, meaning that the bound state has 
dissolved into the continuum. 

With increasing temperature the position of the peak slightly drifts away to the left, meaning that the T(ls) mass 
decreases. The binding energy decreases as well, since the absolute value of the real part of the potential for r, z — >• 00 
increases. In fact, if one subtracts such an asymptotic value from uj, in both figures the peak drifts to the right as 
the temperature increases. However, the bound state does not disappear because the binding energy vanishes, but 
because for T > 400 MeV Landau damping prevents the formation of a bound state. 
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FIG. 4: (Color online) Spectral function of the T(ls) state, divided by —rriQ, at v = for our choice of scales for a s (left panel) 
and for the choice of Refs. [371. [391] (right panel). 



Then, we consider the effect of a non- vanishing velocity. In Fig. [5] we report the spectral functions for the T(ls) at 
T = 250 MeV (left panel) and at T = 400 MeV (right panel) for a few values of the velocity of the plasma. Comparing 
Fig. [3] and Fig. [5] it is apparent that the effect of an increasing velocity - at least qualitatively — is akin to the effect 
of an increasing temperature. Actually, in this case the position of the peak of the spectral function at T — 250 MeV 
does not seem to change at all. But the main effect is that with increasing velocity the height of the peak decreases 
and the corresponding width increases; a behavior that emulates an increase of the temperature of the medium. At 
T = 400 MeV and v = 0, the spectral function has a small peak, which almost disappears at v ~ 0.9. The result at 
this temperature is qualitatively similar to the one observed in [40( at their highest temperature, although one has to 
take into account that we use different reference frames. 

It is interesting to notice that the tendency of the peak to become smaller and wider as the velocity increases, 
changes when going from v = 0.9 to v = 0.99. At v = 0.99 the T(ls) peak is slightly higher and slightly narrower 
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FIG. 5: (Color online) Spectral function of the T(la), divided by -m 2 Q , at T = 250 MeV (left panel) and at T = 400 MeV 
(right panel), for various values of the velocity between the bound state and the plasma. 



than at v = 0.9, as shown for two different temperatures in Fig. [5] This behavior is consistent with the result, 
already discussed in the previous sections, that the bound states become stable at ultra-relativistic velocities. As 
in the particular case considered in the previous section, this behavior can be related to the fact that the effective 
temperature of the plasma is the one given in Eq. (|46p . and therefore for large v the plasma is almost everywhere cold. 
The fact that m£>(z,i>) becomes purely imaginary, implies that the potential ceases to be Yukawa-like and becomes 
oscillatory, as already observed numerically in Ref . [18| • 

Except for this peculiar behavior at v — > 1, both the spectral function analysis and the computation of the width 
through (|32|) show that the width increases as the velocity of the plasma increases, as far as v < 0.9. This is just the 
opposite of the results of Eq. (|29|) in Case I. The reason is that the two results refer to different energy regions, which 
are dominated by different processes. In Case I the thermal width is dominated by gluo-dissociation processes. In the 
present case, the dominant contribution is determined by Landau damping, which is a collisionlcss process. We shall 
further comment on this issue in the Sec. IIV1 



IV. DISCUSSION AND CONCLUSIONS 



In this section we first discuss how the relative velocity v used throughout is related to measurable quantities in 
HIC experiments, like the momentum of the heavy quarkonium state in the lab frame, P^, and the local velocity of 
the QGP, w, in that frame, and make a rough estimate of the importance of the relative motion in the yields. Next we 
compare our result with lattice computations, earlier weak-coupling analysis, and AdS/CFT calculations. We close 
it with the conclusions. 

The clearest experimental signal of the velocity dependence in the in-medium heavy quarkonium properties should 
be in the dilepton yields at fixed rapidity and transverse momentum. In order to have an estimate of the effect, we 
assume that in a central collision the produced medium expands at a constant velocity, w, with respect to the lab 
frame. Typical values for w quoted in the literature are iu|| ~ 1 and w± ~ 0.6 for RHIC and w± ~ 0.66 for LHC. 
We further assume that the system has had enough time to thermalize and that it is isotropic. A heavy quarkonium 
produced with a certain P M in that frame, moves with respect to the plasma with a velocity 



v 



-P°w + ^-w + {P- ^fw) ^l-w 2 

W __V W A ) v 

P 



P°-w 



(51) 



which is the velocity appearing in the formulas of the previous sections. Notice that for a given longitudinal momentum 
this velocity is not totally fixed, it still depends on the transverse momentum and on the modulus of the parallel and 
perpendicular velocities of the plasma in the lab frame, and on the angle ip between w± and P± in the transverse 
plain. The modulus of the velocity can be written as 



{l-w 2 )M 2 



M 2 - 2P°w ■P + (wPy +P 2 



(52) 
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FIG. 6: (Color online) Ratio of the dilepton yields obtained by Eq. (|53[) for the Is state versus P± for typical LHC plasma 
velocities (to ~ mi ~ 0.66). For simplicity we consider P« — 0. 

where M is the heavy quarkonium mass. Assuming a uniform distribution for the angle ip, the modification of the 
dilepton yields can be estimated by the following formula 

Y{v) /' a,r dve -(r(«)-r(«=o))r j (53) 



Y(v = 0) 2tt 



n 



where T(v) is the velocity dependent decay width calculated in the previous sections evaluated for v = v(u>, P) given 
in Eq. (f5"2"j) . r is the lifetime of the QGP, about 10 fm/c for RHIC or about 15 fm/c for LHC, and Y(v) stands 
for the yields. In order to estimate the size of the effect, we display in Fig. [6] the results obtained for Pii = and 
w rv wj_ ~ 0.66 for two different temperatures. We plot the ratio between the velocity dependent yield and the yield 
at v — as a function of P±. The yield has a non-trivial dependence on the transverse momentum, which modifies 
with increasing temperature from a monotonic decreasing behavior at T < 250 MeV to a non-monotonic behavior at 
higher temperature. 

A number of oversimplifications have been done in Eq. (|53p . We have assumed that the heavy quarkonium decays 
in the medium and that the medium temperature and the expanding velocity are constant. These approximations 
should be reasonable if the decay is much shorter than t, and since 1/T(v) ~ 3 fm for T(ls), this seems the case. 
In principle these aspects can be corrected for along the approaches of [41], |42[ or 43 1 (see also (44|)- We have as 
well assumed a constant P for the whole evolution, that is we have neglected the damping of the heavy quarkonium, 
which should be a reasonable approximation because the drift by the expanding medium is expected to be small. We 
have also ignored the velocity dependence in the production mechanism and the contribution of the continuum to the 
yield. Eq. (|53[) is a reasonable approximation if all the above-mentioned corrections factor out in the yield; which 
might not be the case. However, we postpone to future work a more reliable estimate of this quantity. In any case, 
we believe that Eq. (J53"]) , together with Fig. [5] is enough to pinpoint the importance of the relative velocity between 
the heavy quarkonium states and the thermal medium. 

The spectral function of the bottomonium states in a moving thermal bath have been studied with different lattice 
methods in [40j, |45H47| . It is important to take into account that while our computation is performed in the heavy 
quarkonium rest frame, lattice computations are done in the thermal bath rest frame. The spectral function in the 
heavy quarkonium rest frame /OhqC^o) where fco = po — M and p 2 ~ M 2 , is related to the spectral function in the 
plasma rest frame p p i asma by the following equation 

Pplasma(fco) = V{v)pHQ (Wt) > (54) 



where 7 = 1/v 1 — v 2 and 77(11) is a function of the velocity that is not important for the discussion below. If the 
thermal modifications are a perturbation pHQ(fco) in the vicinity of a peak is well approximated by a Breit-Wigner 
distribution 

PHQ(fco) ^(fc„-s/+(F/2)2' (55) 



hence 



Pplasma^o) - ^ _ lE)2 + (7 p /2)2 ■ (56) 
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TABLE I: Comparison of the results reported in th e fig ure 4 of (47J for the r\ c state with the prediction of our Eq. (|42l) . The 
first column p is the momentum in the units used in [47| (~0.5 GeV). The second column is the velocity of the plasma deduced 
by looking at figure 5 in the same reference. The third column is the width obtained by assuming that the spectral function can 
be approximated by a Breit-Wigner distribution and comparing the highest point of the peak with the points where the value is 
half the maximum. Finally the fourth column is the value of the width predicted by using Eq. (|42[) . where, for Tl~ wave (v — 0), 
we have used the corresponding value in the third column (106 MeV). 



This means that even if heavy quarkonium is not modified by the velocity of the plasma in the frame where it is at 
rest, we would still see a modification of the spectral function in the plasma rest frame. This modification will lead 
to an increase of the energy where the peak is located and a broadening of the peak. 

In the case of Sec. HH we can compare with the lattice results of [4fJ. Unfortunately, the velocity range explored in 
[46| was at most of the order of v = 0.2, for which no velocity dependent width change has been observed. This null 
result is compatible with our results taking into account the error in the lattice computations and the fact that their 
temperatures are not very high. Note that the ansatze made in [46[ for the binding energy and the decay width as a 
function of the velocity, based on the hydrogen atom computation in [18| , holds for s-wavc states only according to 
our results for QCD. At zero velocity the lattice results of the same group [48|,[49| turned out to be compatible with 
our results and with those of [13[ . The results of this section for the decay width also appear to be compatible with 
the weak coupling calculation of Ref. [50| at leading order (LO), which also shows a decreasing behavior of the decay 
width with the velocity for small screening masses. At that order, only gluo-dissociation diagrams contribute, like in 
our case. 

In the case of Sec. IIII A[ the results only hold if the thermal corrections can be considered as a perturbation. This 
implies that we can only compare to spectral functions that have approximately a Breit-Wigner form. Because of this 
we cannot compare with the lattice results in [40[ but we can compare with those in [471 ]. By analyzing Figs. 4 and 
5 in [47j we can obtain approximate values of the decay width of the rj c for several momenta. The decay widths we 
obtain from those figures and the corresponding prediction from Eq. (|42[) , which is flavor independent and so holds 
also for charmonium states, are shown in Table |U We observe a similar qualitative behavior, since in both cases the 
width is a non-monotonic function of the velocity, i.e. it increases for low values of v and decreases for larger v, but 
the value of v at which it starts to decrease is lower in [47| than in our case. However, one has to consider that a 
more detailed statistical analysis would be necessary to disentangle possible MEM artifacts from the actual width in 
the plots of Ref. J47J, see [51| 2 . Moreover, the temperatures at which this comparison is done may be too close to 
the deconfinement phase transition for the weak coupling expansion used in our computations to be reliable. The 
equivalent temperature regime for bottomonium (i.e. higher temperature) would be much safer. The results of this 
section agree with the weak coupling estimate of the dissociation temperature in Ref. 52( and also appear to be 
compatible with the contributions to the decay width at next-to-leading order (NLO) displayed in Ref. [5(J. In the 
last reference, it was found that the NLO contribution was much larger than LO one 3 . In our EFT approach this can 
be easily understood if the system is in the kincmatical regime of Sec. IIII1 in which the gluo-dissociation processes 
contributing to their LO are parametrically suppressed. 

A number of analysis on the velocity de pen dence of the screening length have been carried out for strong coupling 
using the AdS/CFT approach [53l - l59l ] (see (6Q| for a review). It is not straightforward to compare these results to ours, 
as they do not obtain an imaginary part in the potential. This would be a first important difference. Furthermore, 
in momentum space, what plays the role of the screening mass for us is the complex, angle and velocity dependent, 
Debye mass m£i(v,0), see Fig|3] , which translates into a non-trivial potential in coordinate space for which no simple 
analytical form has been found. Hence we cannot make further statements on this respect than those already made in 
Ref. [18| 4 . However, we can certainly compare with the two AdS/CFT calculations of the heavy quarkonium spectral 
function at non- vanishing velocity we are aware of [6ll . |62[ . These spectral functions qualitatively agree with ours 



2 We thank Masayuki Asakawa for pointing this out to us. 

3 There are a number of approximations in the NLO calculation of [50 
discussion. 

4 In formula (92) of 



lq . m,£)(v,8) should read \m,£)(v,8)\. 



in particular Pauli blocking is ignored, see [221 ] for a recent 
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in Case II at moderate velocities, in the sense that the bound state peaks become smaller and wider as the velocity 
increases. Let us finally remark that we observe in the ultrarelativistic limit an oscillatory behavior of the potential 
rather than an exponential damping, that would lead to the stabilization of the bound states, which is not observed 
in the AdS/CFT approach. 

In summary, we have analyzed heavy quarkonium states moving in a weakly coupled QCD plasma. In the Case I, 
corresponding to the hierarchy wiq 3> 1/r 3> T 3> E 3> tod, we have found that the thermal decay width decreases as 



the velocity increases, like in QED [18j]. The decay width is in this case dominated by gluo-dissociation processes [21 
However, unlike in QED, the thermal energy shift becomes velocity dependent, except for the s-wave states. In the 
Case II, corresponding to the hierarchy ttiq ^$> T ^$> 1/r , mu 3> E, we have found a different behavior for the decay 
width, namely it increases as the velocity increases, except for ultra-relativistic velocities for which it starts decreasing 
again. This non-trivial behavior was overlooked in Ref. [18| . The decay width is in this case dominated by the Landau 
damping. Putting all together, we conclude that the decay width depends in a nontrivial way on the temperature and 
on the velocity, which complicates the interpretation of HIC experimental data, as we tried to illustrate by Fig.jGl Our 
results are consistent and in qualitative and semi-quantitative agreement with the few available lattice data j46l |47| , 
and also appear to be compatible with the weak coupling analysis of refs. |50ll52|. 
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Appendix A: General Framework 

In this appendix we briefly review the general framework used to take into account the effect of a moving thermal 
medium. A more detailed discussion can be found in [18[ . 

We shall assume that the plasma (or black-body radiation) is in thermal equilibrium at a temperature T. Since we 
are considering the reference frame in which the plasma is moving with a velocity v, the particle distribution functions 
are given by 

If.b^K) = — L— - , (Al) 

elf v ± i 



where the plus (minus) sign refers to fermions (bosons). In the reference frame where the thermal bath is at rest 



fi^k^ = &, while in a frame where the plasma moves with a velocity v we have that 



r 5" = J(l ) t») = Y. ( A2 ) 



where 7 = l/\/l — v 2 is the Lorentz factor; the latter frame has been successfully used in the past, for example in 29f. 
Studying a bound state in a moving thermal bath is akin to study a bound state in non-equilibrium field theory 63J; 
in that case the Bose-Einstein or Fermi-Dirac distribution functions are substituted by a general distribution, which in 
our case will be the boosted Bose-Einstein or Fermi-Dirac distribution functions reported in Eq. (|A1|) . For a thermal 
medium formed of massless particles, taking into account that in non-equilibrium field theory the collective behavior 
always enters through on-shell particles or antiparticles, we have (in the case of particles) that 

, 1 — vcosO 
^k u = k ; , A3) 

where k — \k\ and 6 is the angle between fc and v. The distribution functions in Eq. (jAlj) can now be written as 

f F Mk,T,e,v) = ek/T J e>v)±1 , (A4) 

where we have defined the effective temperature 



TVl - v 2 

Teff(&,v) = ^ V -, (A5) 

1 — v cos 9 
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which is plotted in Fig. [7]for few values of v. Fig. [7]hclps to clear away the misconception that a bound state moving 
with non- vanishing velocity in a thermal bath feels a higher temperature. Indeed, the effective temperature is in most 
of the directions smaller than T; for v ~ 1 we find that T e ff(9, v) > T only for < 6 < \/2(l — w 2 ) 1 / 4 . Intuitively, the 
dependence of the effective temperature on v and 9 can be understood as a Doppler effect. 



\ 



v=0.1 

v=0.5 

v=0.9 

v=0.99 



■J ? 2 - \ V 



ji/4 



ji/2 

e 



3ji/4 



FIG. 7: (Color online) Effective temperature divided by T, as a function of the angle between the vectors k and v, for 
v = 0.1, 0.5, 0.9, 0.99. 



While at v = it is clear that the thermal medium introduces a new scale T in the problem, it is not clear a priori 
how many scales a moving thermal nrcdium introduces. This can be understood by using light-cone coordinates. We 
choose v in the z direction and define 



fc+ = fcn 



Then, we have that 



where 



to 



and 



k- = fc n - fa 



^ = U k -± + k - 
p 2 I T T _ 



(A6) 



(A7) 



T + =T 



1 + v 



1 



and 



T_=T 



1-v 



1 



(A8) 



Therefore, in light-cone coordinates, it becomes explicit that the distribution function actually depends on two scales, 
T + and T_, Obviously, for any value of v, one has that T + > T > T_, and moreover T + corresponds to the highest 
temperature measurable by the observer, while T_ corresponds to the lowest temperature measurable by the observer. 
In this work we consider always that T + and T_ arc of the same order of magnitude. Even though this is not so for 



very large velocities, in all the cases considered in 
correct even for v — > 1. 



18j we found the results obtained assuming T + ~ T_ were indeed 



Appendix B: Numerical checks 

The operator-splitting is a powerful method for solving partial differential equations |64| . The idea is to separate 
a complex differential equation in various simpler equations and to solve them with a discretization method. For 
illustrative purposes we compare the results obtained for a hydrogen-like atom with two different potentials, namely 
the Yukawa potential 



Vy(x) 



a 



\J x 1 + y 2 + z 2 



-tod y/x 2 +y 2 +z 2 



(Bl) 
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and the potential reported in [18| for vanishing velocity of the thermal medium. These analyzes could be done 
employing a numerical code with spherical symmetry, however using cylindrical coordinates allows us to check the 
numerical procedure employed in the general case of non-vanishing velocity. Indeed, the present discussion can be 
generalized to any potential with cylindrical symmetry. 

The Schrodinger equation in cylindrical coordinates, (z,r,ip), for s-wave states is given by (for simplicity we set 
m e = 1 , a = 1 ) 

du(t,r,z) fid 2 Id 2 Id 1 Tr , A , 

being u(t,r, z) — rip(t,r,z); in order to compute the spectral function we use the initial condition u(0,r, z) = 
—rS(r)S(z) and boundary condition uit, 0, 0) = 0. We separate the Hamiltonian in one term containing the po- 
tential and derivatives with respect to z and a second term with only derivatives with respect to r: 

(B3) 

(B4) 

and we solve the corresponding Schrodinger equations numerically using the Crank-Nicolson method [65j, meaning 
that the equations are discretized as follows, 

l+ l -H iat \u n+1 / 2 = (l - % -H ia \ u n , (B5) 

1 + iffeat) u n+1 = (l - l -H 2 a t ^j u n +^ 2 , (B6) 

where a t is the temporal lattice spacing and n indicates the discretized time step. The equations are solved recursively: 
At the n-th step the wave-function evolves from time nat to (n + l/2)ot according to Eq. (|B5| . In the next step 
from (n + l/2)at to (n + l)dt according to Eq. (|B6[) . In Eqs. ()B5|) and (|B6|) also the space coordinates are discretized 
z = la s and r = ja s , where a s is the spatial lattice spacing and I, j are integers. 
The discretized initial condition reads 





1 d 2 
-55. + V(r,.), 


2 — 


Id 2 Id 1 

_ 2cV2 + 27cV ~ 2^2 



rS 2 (r)6(z) = ^ J d 2 p r e l ^ r J°° d Pz e^ z 



dp r — Ji (p r a s ) J (j'p r a s ) / dp z e tPz la ° 



4?H J a a s J-K/a, 

j r 

= ? s io du Ji(u) J (j u) , (B7) 

where the Bessel functions of the first kind, J a , are used instead of the trigonometric function of [371 ] to improve the 
convergence. Once the discretized wave functions arc obtained, the corresponding spectral function can be calculated 
by means of Eq. (JT7)) . 

We consider first the Yukawa potential. The binding energy of the various states is known with great accuracy, see 
e.g. [66M68I ]. and it is also known that sequential dissociation of the bound states takes place with increasing values 
of A = mjjao, where a® is the Bohr radius. On the left panel of Fig. [5] we report the spectral functions obtained 
with the splitting method for various values of A; the vertical lines correspond to the numerical values of the binding 
energy obtained in [68[ for the Is state. As expected, at mp = (green line) we find a peak at u> ~ — 0.5m e a 2 , 
corresponding to the Is state of the standard hydrogen atom. The 2s state at cj ~ — 0.125m e a 2 is also visible, but 
the corresponding height is suppressed. In principle, at any bound state should correspond a Dirac delta-function, 
however the discretization procedure can only lead to a finite peak. The height of the peak is proportional to the field 
strength, that is to the strength of the corresponding interaction channel, which explains the fact that the 2s peak is 
smaller than the Is peak. 

With increasing values of the Debye mass the peak of the Is state moves to higher values of energy, meaning that 
the corresponding binding energy decreases. Although the peak height decreases, for the same reason explained above, 
note that no appreciable broadening of the spectral function appears, meaning that the numerical procedure does not 
produce a fictitious increase of the width. Indeed, in this simple model (and in any model with a real potential) the 
dissociation happens when the peak of the spectral function approaches zero. Our numerical results indicate that the 
Is state of the Yukawa potential dissociates at A ~ 1.2, in good agreement with the numerical results of [67ll68l|. 
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FIG. 8: (Color online) Left panel: Spectral functions obtained for the hydrogen-like atom with the Yukawa potential (I 
for various values of A = rnD/(m e a)- The numerical algorithm slightly overestimates the binding energies obtained in ( 
corresponding to the vertical dotted lines. Right panel: Spectral functions obtained with the potential reported in [la 
vanishing velocity. 



Regarding the 2s state, at A = 0.1 (blue line) it is still visible, with binding energy u> ~ — 0.05m e a 2 (in good 
agreement with the results of [63, [68|), but for larger values of A the 2s state is no more visible, although it is known 
that it only dissociates at A ~ 0.31, see [63, |68[. The reason, as explained above, is that the corresponding peak is 
very small and cannot be identified with the used numerical accuracy. 

In summary, from the analysis of the Yukawa potential, we conclude that the algorithm correctly reproduces the 
binding energy of the Is state at any value of the Debye mass, but the decrease of the peak height observed in the left 
panel of Fig. [5]for increasing values of ttid is an artifact due to the combined effect of the numerical discretization and 
of the reduction of the strength of the interaction channel. Remarkably, the algorithm does not produce a fictitious 
width. The analysis of the dissociation of the excited states for this model with the present method is problematic, 
because of the reduction of the peak height. 

On the right panel of Fig. [51 we report the spectral functions obtained with the potential reported in [18[ considering 
vanishing velocity. This potential has an imaginary component for any non-vanishing value of the temperature. At 
T = 0, i.e. for A = 0, the potential is real and Coulombic, and the standard peaks of the hydrogen atom for the Is 
state at ui ~ — 0.5m e a 2 , the 2s state at w ~ — 0.125m e cv 2 and the 3s state at w ~ — 0.05m e a 2 are reproduced with 
a good accuracy. As before, the peaks of states with high principal quantum number are suppressed. Increasing the 
temperature, the binding energy of the Is state decreases, but the corresponding spectral function not only moves to 
higher energies, it also becomes wider. From the insight gained in the analysis of the Yukawa potential, we conclude 
that the broadening of the peak is due to the imaginary part of the potential and not to the numerical procedure. 
Moreover, improving the discretization procedure, we checked that the reduction of the peak height of the Is state 
(at nonvanishing temperature) is not an artifact, but it is instead a genuine effect related to the imaginary part of 
the potential. The reason is that with a finite imaginary potential, the spectral functions are not delta-functions, but 
smoother functions which can be resolved with the used discretization method. 
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